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The counting rule for Nambu-Goldstone modes is discussed using Mori's projection operator 
method in nonrelativistic systems at zero and finite temperatures. We show that the number of 
Nambu-Goldstone modes is equal to the number of broken charges, Q a , minus half the rank of the 
expectation value of [Q a , Qb\- 
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Introduction — Symmetry and its spontaneous break- 
ing are of basic importance for understanding the low- 
energy physics in many-body systems. When a contin- 
uum symmetry is spontaneously broken, there exist zero 
modes called Nambu-Goldstone (NG) modes [1-3], e.g., 
pions in hadron physics, spin waves in a ferromagnet, 
or phonons in a crystal. However, the number of NG 
modes associated with the spontaneous breaking appear 
in a nonrelativistic system is an unsolved problem [4] . In 
the case of relativistic systems (more precisely, Lorentz 
invariant systems), the number of NG modes coincides 
with the number of broken symmetries, and their dis- 
persion is linear, i.e., u> = |fe|, where w and k are the 
energy and momentum of the NG modes. In contrast, in 
the case of nonrelativistic systems, the number of NG 
modes is not necessarily equal to the number of bro- 
ken symmetries, and the dispersion may be nonlinear. 
Nielsen and Chadha (NC) [5] classified NG modes into 
two types: type-I and type-II NG modes whose ener- 
gies are odd and even powers of \k\, respectively. They 
showed that the number of type-I NG modes plus twice 
the number of type-II NG modes is greater than or equal 
to the number of broken symmetries. In known realistic 
examples of type-H NG modes such as spin waves and 
NG modes in the Kaon condensed phase [6, 7], the in- 
equality is saturated. The presence of an expectation 
value for charge densities implies a type-H dispersion for 
NG modes, as was clarified by Leutwyler in an effective 
field theory [8]. Schafer et al. pointed out that the ex- 
pectation value of [Q a ,Qb] rather than charge densities 
themselves plays important role, and if they all vanish, 
the number of NG modes is equal to the number of bro- 
ken symmetries, where Q a are the broken charges [7]. 
Recently, this theorem was generalized to the inequality, 



N B s-N NG < -rank([Q a ,Q fc ]), 



(1) 



by Watanabe and Brauner (WB) [9]. They conjectured 
that the inequality in Eq. (1) is saturated for general sys- 
tems. These results were obtained at zero temperature. 
In this Letter, we generalize the Nambu-Goldstone the- 
orem in relativistic systems to that in nonrelativistic sys- 
tems at zero and finite temperatures, and we prove that 
the inequality in Eq. (1) is saturated. We give an al- 



ternative definition of type-I and type-II NG modes by 
using a matrix ([Q a ,Qb]), arid show that the number of 
type-H NG modes is equal to rank([Q a ,Qf,])/2, and the 
NC inequality is also saturated. 

For this purpose we apply Mori's projection opera- 
tor method [10] to systems in which spontaneously sym- 
metry breaking occurs. In this approach, a Hamilton 
(Langevin) equation for NG modes is derived at zero (fi- 
nite) temperature, in which the expectation value of the 
commutation relations for NG fields and charge densi- 
ties become the Poisson brackets for their fields, i.e., NG 
fields and charge densities are canonical variables. The 
Hamiltonian formalism for NG modes in nonrelativistic 
systems at zero temperature was discussed in an analy- 
sis of specific models [11]. We develop the analysis in a 
model independent and nonperturbative way. 

In the following, we consider the case that the trans- 
lation symmetry is not spontaneously broken, and that 
the broken charges do not explicitly depend on space- 
time variables, i.e., the conserved charge densities are 
uniform [9], n a (t,x) — e lPx n a (0, 0)e~ , where P M are 
the spacetime-translation operators. 

Mori's projection operator method — Here, we briefly 
review Mori's projection operator method [10], which is 
a powerful tool for discussing the low-energy excitations, 
and essential for proving our theorem for nonrelativis- 
tic systems. For readers who are unfamiliar with this 
method, see, e.g., Ref. [12-16] for more details of this 
method. 

Consider a set of operators {A n (t,x)} corresponding 
to the soft modes, chosen to be NG fields and broken 
charge densities. First, we define the thermal average 
of an arbitrary operator as (0) = tr p cq O, where we 
choose p C q as the canonical-ensemble-density operator, 
p cq = cxp(— /3iJ)/trexp(— (3H), with the Hamiltonian H 
and the inverse temperature f3 = 1/T. One may consider 
the grand canonical ensemble, by replacing H with H — 
fxN, where /i and N are the chemical potential and the 
conserved charge operator, respectively. 

Next, we introduce an inner product for arbitrary op- 
erators C>i and 02: 



(0 x ,0 2 ) = i [ dT{e TH O x e- TH 0\), 

P Jo 



(2) 



which satisfies positive definiteness, (0i,0i) > 0, and 
conjugate symmetry, (0i,02) = (C2,Ci)*- Using this 
inner product, we define a metric: 



g nm (x-y) = (A n (0,x),A m (0,y)). 



(3) 



We also define g (y — z) as the inverse of g nm ( x — y), 
i.e., 

J d 3 yg nm (x - y)g ml (y - z) = 5 n l d^(x - y). (4) 

Here, we used Einstein notation: if an index appears 
twice in a single term, once as a superscript and once as 
subscript, a summation is assumed over all of its possible 
values. An operator with an upper index is defined as 



A n (t,x)= d i yg nm {x-y)A m {t,y). 



(5) 



Since the metric is the two-point function of A n , the in- 
verse metric coincides with the second derivative of the 
effective action, (3r(A n ), with respect to A n : 



g ml (x-y) = 



s 2 pr(A n ) 

5Ai(y)dAL(x)' 



(6) 



where the effective action is given by the Legendre trans- 
formation of the generating functional W(J n ): 



r{A n ) = W(J n ) - / d 6 xJ m {x) 



t 5W{J n ) 
8J m (x) ' 



with 

-/W(J") 



trexp 



0H + / d 3 xA n (0,x)J n (x) 



(7) 



(8) 



In order to decompose the degrees of freedom into the 
slow modes and others, we introduce the projection op- 
erator V acting on a field B, given by 



VB(t, x) = I d a yA n (0, x)(B(t, y),A n (0, x)). (9) 

We also define Q = l-V. These satisfy V 2 = V, Q 2 = Q, 
and QP = VQ = 0. 

Now, let us derive a generalize Langevin equation from 
Mori's projection operator method by reconstructing the 
Liouville equation, doA n (t, x) = i£A n (t, x), where CO = 
[H, O] . The formal solution of the Liouville equation is 
obtained as A n (t,x) — e lCt A n (0,x). 

First, we decompose <9 e into 



d e l " = e'"i£ = e tLt ViL + e lLt QiL 



iCtr, 



ACt t 



(10) 



Next, consider the Laplace transform of exp(i£i), which 
can be decomposed into 



1 



1 



1 



z — iC z — QiC z — iC 



ViC 



1 



z- QiC 



(11) 



The inverse Laplace transform leads 

e iCt = e QiCt + f dse iC{ - t - s) PiCe QiCs . (12) 

Substituting Eq. (12) into the second term in the right 
hand side of Eq. (10), we obtain the operator identity 

d e tCt = e tCt ViC+ f dse lC ^-^ViCe QlCs QiC+e QlCt QiC. 
Jo 

(13) 

Acting Eq. (13) to A n (0, x), we obtain the generalized 

Langevin equation of motion [10], 

d t A n (t,x) = / d 3 yi!2 n m (x ~ y)A m (t,y) 

- ds d 3 yK n m {t -s,x- y)A m {s, y) 



+ R n (t,x). 



(14) 



Here if2 n m (x — y) is the streaming term, K n m (t—s, x—y) 
is the memory function, and R n (t,x) is the noise term 
defined as 

iQ n m {x -y) = (iCA n (0, x), A m (0, y)) 



P 



([A n (0,x),A^(0,y)]}, (15) 



K n m (t-s,x-y) = -e(t-s)(iCR n (t,x),A m (s,y)), 

(16) 

R n (t, x) = e ltQC QiCA n (0, x), (17) 

where 0(t) is the step function. Note that this equation 
is an operator identity and thus Eq. (14) is equivalent to 
the Liouville equation. The noise term is orthogonal to 
A m (0,x), i.e., (A m (0,x),Rn(t,y)) = 0, so that it drops 
out when the two point function (A n (t, x), A m (Q, y)) is 
considered. Therefore, the dispersion relation for A n can 
be determined by the linear equation with if2 n m (x — y) 
and K n m (t,x — y). In momentum space, we can write 
the equation as 



d t A(t,k) =iO(k)A(t,k) 



/ dsK(t - s, k)A(s, k) + R(t, fe), 



(18) 



in matrix notation. This is the basic equation used in 
the our analysis. 

Nambu-Goldstone theorem — Let us start with the 
Nambu-Goldstone theorem for the effective potential [3] . 
Consider a set of conserved charges, 



Q a = d 3 xn a (t,x), 



(19) 



which commutes with the Hamiltonian: [Q a ,H] = 0. 
When the symmetry is spontaneously broken, there ex- 
ists a matrix, 



i([Q ai ^i(t,x)}) = [M n ^} c 



(20) 



such that M n< p is regular, where i and a run from 1 to 
the number of broken symmetries TVbs • We assume that 
(f>i(x) are not conserved charges. Without loss of gener- 
ality, 4>i can be chosen to be a real field. 

Here, let us consider the effective potential defined as 
V(#i,JV„) = r{A n )/V by choosing A n = (*,-,#„) that 
contain NG fields and broken charge densities, where V 
denotes the space volume. We parametrize <P = (ipiiVj) 
and N = (n a ,n' b ), where <f>i (ipj) and n a (n' b ) are (non-) 
NG fields and (un-) broken charge densities. 

The effective potential is invariant under the infinites- 
imal symmetry transformations: 

8$i = -ie a {[Q a ,$i}) = e a [M Nf ] ai , (21) 

SN b = -ie a ([Q a ,N b }) = e a [M NN ] ab , (22) 

that is, 

SV SV 

[M N 4,] at — + [M NN ] ac j^ = 0. (23) 

This is an identity of the effective potential. If <Pi and 
N a belong to linear representations of a Lie group, they 
become [M N<P ] ai = -i[T a &]i, and [M NN ] ac = f a b C N c , 
where T a and f a b C are the generator and the structure 
constant, respectively. Differentiating Eq. (23) with re- 
spect to <Pi or N a , we obtain 



5[M, 



N&\ai 



SV 



8$, 5$i 



[M, 



N&\ai 



+ [M NN ] ac 



5[M. 



NN\ 



SV 



SN h SN n 



[M 



N&\ 



5 2 V 
S$iS$j 

S 2 V 
5N C 5$ 3 

6 2 V 
'SQiSNt 

S 2 V 



= 0, 



+ [MNN] -lNJN b =°- 



(24) 



(25) 



At the stationary point, the first terms in Eqs. (24) and 
(25) drop, and only terms with NG fields and broken 
charges survive. Equations (24) and (25) become 



M n4> V^ + M nn V 



n<j> 



o. 



M n(/> V* n + M nn V nn = 0, 
(26) 
where we used matrix notation, and V a/3 = S 2 V / (SaS f3) 
with a, — {4>i,n a ). Therefore, one finds that 

(4>\\nh,') = ([M^n]ib, [M nn ] ab ) are eigenvectors of V af} 
with the vanishing eigenvalue. The number of eigen- 
vectors is equal to the number of broken symmetries. 
In Lorcntz invariant systems, M nn — because non 
Lorentz scalar operators cannot condense. In this case, 
the inverse of propagator at fc M — 0, V^, vanishes from 
Eq. (26), which implies the propagator of <j> have poles 
at k 2 = 0. The number of NG modes coincides with the 
number of eigenvectors, i.e., the number of broken sym- 
metries, and their dispersion is ui — \k\. This is the NG 



theorem in Lorentz invariant systems. In the case of non- 
relativistic systems, however, the number of NG modes 
may differ from the number of broken symmetries. 

To derive the counting rule for NG modes in nonrela- 
tivistic systems, we choose the operators in Eq. (14) as 
A n (x) = (4>i(x),h a (x)) with ^(x) = <f>i{x)- (<f>i(x)) and 
h a (x) = n a (x) — (n a (x)). We assume that no other zero 
modes that couple to NG modes exist in the system. If 
this is not the case, we must add a field coupled to the 
zero mode to A n as a slow variable and solve the coupled 
equations. 

Now, let us consider the Laplace transformation of the 
memory function, K(z, fc), which gives the dissipation of 
NG modes. We are interested in the behavior of the mem- 
ory function at low energy and low momentum, so we 
take z — > and fc — >• 0. Since Cn a (z,k) = — fc ■ j a (z,k), 
where j a (z,k) is the current operator, K nn (z,k) ~ fc 2 , 
K n< p(z,k) ~ fc, they vanish at fc = 0. Only K,f,^(z,k) 
can survive at fc = 0. K c j >( j > (z, 0) can be expanded as 

0K H (z, 0) = SM H + L H + z5Z H + 0(z 2 ). (27) 

The first and second terms SM^ and L^, is called the 
Onsager coefficients. Both are real, and SMT, = —SM$$ 
and LT, — iX, are satisfied. SM^ gives a correction 
of M^0. Lfij,, which represents dissipation of the field 
is related to a spectral function at the zero frequency, 
and it vanishes at zero temperature because there is no 
spectrum at the zero frequency and the zero momen- 
tum. Since M^(k) and K^(z,k) appear as a com- 
bination Mss(k) — K^(z,k) in the Langevin equation 
[see Eq. (18)], SM^ can be renormalized to M^ by the 
shift M 00 -► M 00 - SMm and 0K^ -t 0K^ - 5M^. 
The third term SZ^ denotes the correction of the wave 
function. In the following we drop SZ^, which will be 
justified later in the analysis of zero modes. First let us 
consider the case at zero temperature dropping the mem- 
ory function. Then, the equation of motion becomes 



d 



4>(t,k)\ _ fiQ/(k) i^ n (k) 



n(t, fc) 



iHJik) if2 n n {k) 



4>(t,k) 

h(t, fc) 



(28) 



with iHj(k) = M ay (k)r^(k), where M a/3 (fc) 
and r al3 (k) is defined as the Fourier transform of 
-i([a(0,x), 0(0,0)]) and 5 2 r/(Sa(x)S0(O)), which coin- 
cides with M a p and V"' 3 at fc = 0, respectively. Equa- 
tion (28) has a form of a Hamilton equation. M a p(k) 
can be identified as the Poisson bracket for a and 0. We 
assume &etM a p{k) ^ 0, which ensures 4> and n are in- 
dependent canonical variables. r Q/3 (fc) corresponds to 
the hessian of the Hamiltonian, whose eigenvalues are all 
positive at fc ^ 0, and TVbs of them becomes zero at 
fc = 0, i.e., rank(V"^) = A^bs- The eigenmodes are ob- 
tained by det(iuj + if}) = 0. u have always real ±E n pairs 
(E n > 0, n = 1, • • • , Nbs) since this is the Hamiltonian 
system with nonnegative hessian. 



At fc = 0, i74 (O) and iQ n n {0) vanish from Eq. (26), 
which is equivalent to the conservation law dQ a /dt = 0. 
We also obtain 



iJV(O) = F~\ tf2/(0) = -F- X G, 



(29) 



where F" 1 = (M$ n - M^M^)- 1 M nn )V nn can be 
identified as the inverse decay constant (matrix) for the 
NG modes up to a renormalization factor, and G — 
M nn M^ n ~ x . F~ x is a regular matrix, which can be con- 
firmed as follows: Consider V = M a sV Sl M 1 p. Since 
M a p is a regular matrix, rank(V') = rank(V"' 3 ) = Abs 
is satisfied. V' is explicitly obtained as 



V 



F-^-GM^ + M^) 




(30) 



Then, rank(V') = rank(F- 1 (-GM 00 + M n0 )) = A BS ; 
therefore, rank(F _1 ) = Nbs- The rank is equal to the 
matrix size, so that -F -1 is a regular matrix. 
The equations of motion become 



d j>(t,0) = -F- 1 GJ>(t 1 0) 
d n{t, 0) = 0. 



F^nftO), 



(31) 

(32) 



Then, the eigenvalue equation is det(iw + ifi) — 
(iuj) NBS det(«w - F- X G) = 0. The number of 
massive modes coincides with rank(F~ 1 G)/2 = 
rank(F- 1 M„ n M ( ^ n 1 )/2 = rank(M„„)/2. Therefore, we 
arrive at 



Nbs ~ A NG = - rank(M„„). 



(33) 



Equation (33) is nothing but the WB conjecture. 

At a finite temperature, the Onsager coefficient, L^, 
contributes to the equation of motion. This effect can 
be included by replacing M^ with M^ — L^. In this 
case, the above counting rule does not change as long as 
det M a p ^ 0. 

Next, let us classify the NG modes. Differentiating 
Eq. (31) by d , we obtain 



d^{t,0) = -F- l Gdo$(t,0). 



(34) 



We define the NG fields satisfying G<j>i — as type-I NG 
fields, and others that are linearly independent of cf>i as 
type-H NG fields 0n ■ Obviously, the number of type-I 
NG fields is equal to Nbs — rank(M„„), which coincides 
with the number of type-I NG bosons, N 



typc-Ij 



because 



<9o</>i = 0. On the other hand, the number of type-H 



NG modes is equal to rank(M n „)/2 
A NG = At y pc-i + A typo _n , we find 



A, 



typc-II 



Abs = Atype-i + 2A typo _n , 



Since 



(35) 



from Eq. (33) . This is the saturation of the NC inequal- 
ity [5]. Equations (33) and (35) are the main results 



in this Letter. Here, let us discuss whether the correc- 
tion of the wave function SZ^ in Eq. (27) is relevant 
to the counting rule of the zero modes. bZ^ changes 
d$ to (1 - <5Z 00 G T V""G)<9 2 . For the type-I NG modes, 
this correction is irreverent because G<t>\ — 0. For the 
type-H NG modes, the second derivative is higher order; 
therefore, SZ^ is irreverent to the counting rule of NG 
modes. 

Explicit breaking and mass formulae — Here we discuss 
the masses of the NG modes when an explicit breaking 
term is added. Suppose the symmetry is explicitly broken 
by an interaction term 5V — <Pih l . In this case, the 
contribution from the explicit breaking term is obtained 
from Eq. (24): 



^Jm 



6$i 



h J , 



(36) 



where we used SV/S'Pi = h % . Therefore, the equation of 
motion becomes 

d 2 <Xi, 0) = -F-i-Gdofo, 0) + F-HOJiOMt, 0). (37) 

For general cases, it is not easy to derive the disper- 
sion relation from Eq. (37), so here we consider a spe- 
cial case: the symmetry breaking term does not mix the 
type-I and the type-IE fields, so we can decompose the ex- 
plicit breaking term to F _1 z/?„ (O) = [F~ 1 if2 n % ypc .i ® 
[F~ 1 if2 n 'P] t y pB _]i . Then, the equation of motion for <pi 
reads 



a o 2 0i(t,o) = [F-HQjjtype-jM^o), 



(38) 



and thus, the squared mass matrix for type-I NG modes 
becomes w 2 ypc _j = — [F^ l iQ n ^] typc .\ = 0(h), which cor- 
responds to a generalized Gell-mann-Oakes-Renner re- 
lation [17]; m 2 = — TTiqfy^xJ j where (ij)ip), xJs > an d 
m q correspond to M,p n , V"™, and i£2„ : . , respectively. 
On the other hand, the equation of motion for 0n is 

F-iGdoifa (t, 0) = i[F- l irt n *] 1y p-n&i (t, 0), (39) 

where we dropped the cJg^n (t, 0) term because it is 
higher order in h. Since the do ~ h in Eq. (39), the 
squared mass becomes mjL n = 0(h 2 ). 

The NG modes of type-I and type-H have different type 
of mass formula: the squared mass is proportional to the 
h and h 2 for type-I and type-II NG modes, respectively. 
The order counting of the mass matrices for type-I and 
type-H NG modes does not change even if the mixing 
term between their modes exists in F~ x iQ n ® [18]. For 
a typical case at finite momentum, the order of the mo- 
mentum dependence corresponds to h ~ fc 2 , then the 
dispersion of type-I NG modes is u) ~ |fc|, and that of 
type-H NG modes is u ~ fc 2 . In this case, the classifi- 
cation of type-I and type-H in this Letter is the same as 
that by Nielsen-Chadha. 



Summary and outlook — We have derived the counting 
rule for the NG modes, Eqs. (33) and (35), in nonrel- 
ativistic systems at zero and finite temperatures. Our 
method is model independent; the details of the the- 
ory are reflected by the expectation value of commuta- 
tion relations for the NG fields and the broken charges, 
as well as the second derivative of the effective action. 
To derive the counting rule, we employed the following 
assumptions: (1) The broken charges are uniform; (2) 
Translation symmetry is not broken; (3) det M ni p ^ 
and det M a p ^ exist. The assumption (1) implies that 
the local operator does not explicitly depends on space- 
time variables. The assumption (2) is necessary that NG 
modes are to be eigenmodes of momentum. The assump- 
tion (3) ensures that <f> and n are independent canonical 
variables. In this Letter we assumed that det([Q a , <f>i\) ^ 
and 4>i are not conserved charges. If a set of conserved 
charges Q' a exist such that ([Q' a , 4>i[) = ([Q' a ,Qb]) = 
for all <f> t and Q b & Q' a ) but ([Q' a) Q' b ]) + 0, Q' a are 
also type-II NG fields [19]. This can be derived using the 
same reasoning used in this work. The number of these 
modes is equal to ra,nk{Q' a ,Q' b ]/2. Therefore, our results 
in Eqs. (33) and (35) remain unchanged. 

Our method can also apply to the systems where space- 
time symmetry is spontaneously broken; this is beyond 
the scope of this work. 
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Yazaki for useful discussions. This work was supported 
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Note added: After finishing this work we became aware 
of similar work by Watanabe and Murayama [20]. Al- 
though their analysis is different from ours, obtained re- 
sults are consistent with our work. 
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